Propagation of entangled light pulses through dispersing and absorbing channels 
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The problem of decorrelation of entangled (squeezed- vacuum-type) light pulses of arbitrary shape 
passing through dispersive and absorbing four-port devices of arbitrary frequency response is studied, 
applying recently obtained results on quantum state transformation [Phys. Rev. A 59, 4716 (1999)]. 
The fidelity and indices of (quantum) correlation based on the von Neumann entropy are calculated, 
with special emphasis on the dependence on the mean photon number, the pulse shape, and the 
frequency response of the devices. In particular, it is shown that the quantum correlations can 
decay very rapidly due to dispersion and absorption, and the degree of degradation intensifies with 
. increasing mean photon number. 
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00 ! I. INTRODUCTION 

£f*l • Quantum-state entanglement of spatially separated systems is one of the most exciting features of quantum me- 
| chanics jjj. In particular in the rapidly developing field of quantum information processing (quantum teleportation 
P), quantum cryptography 0, and quantum computing Q), entanglement has been a subject of intense studies. 
Recently, continuous- variable systems have been of increasing interest p|-|l6[]. 

Since entanglement is a highly nonclassical property, it is expected to respond very sensitively to environment 
influences, and thus it can decrease very fast. The mechanisms of decoherence are worth to be studied in detail 
]l7| . p^| , because they delimit possible applications. In particular, optical pulses prepared in entangled states typically 
propagate through optical fibers and/or pass optical instruments, such as beam splitters, mirrors, and interferometers. 
All these devices are built up by (dielectric) matter that always gives rise to some dispersion and absorption. As a 
result, the initially prepared quantum coherence is destroyed, and the question arises of what is the characteristic 
scale of quantum decorrelation. 

In this paper the quantum decorrelation of two initially entangled light pulses that pass through optical devices is 
studied. The pulses are regarded as being nonmonochromatic modes of arbitrary shape, and the devices are regarded 
as being dispersing and absorbing four-port devices of arbitrary frequency response. The underlying theory of quantum 
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O^' state transformation has been developed recently |19|. It is based on a quantization procedure for the electromagnetic 
field in dispersing and absorbing inhomogeneous dielectrics [po|-p4|, which is consistent with both the dissipation- 
fluctuation theorem and the QED canonical (equal-time) commutation relations. One interesting aspect of the theory 
is that given the complex refractive-index profiles p5| of the devices, which can be determined experimentally, the 
parameters relevant to the quantum-state transformation can be calculated without further assumptions. 

In Section || the basic relations for describing the pulse propagation and the associated quantum-state transforma- 
tion are given for the case when the pulses are initially prepared in a two-mode squeezed vacuum state. The fidelity 



and some characteristic measures of (quantum) correlations of the pulses are calculated and discussed in Section III 



and some concluding remarks are given in Section IV 



II. PULSE PROPAGATION AND QUANTUM-STATE TRANSFORMATION 

Let us consider the propagation of two nonclassical light pulses through lossy four-port devices of given complex 
refractive-index profiles P5J. Regarding the pulses as nonmonochromatic modes, we may define annihilation operators 
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for the two pulses in terms of the annihilation operators associated with the monochromatic modes that form the 
pulses Mftf], 



/>oc 

a[r)] = / dui rf (u>)a(ui), 
Jo 

poo 

d\fj\ = I duj fj* (ui)d(ui), 
Jo 

where 7](u>) and fj(u>) are normalized functions that describe the pulse profiles, 

r- OO 

\\r,f = (r,\rj) = / dojri*(u,)r ] (uj) = l, 
Jo 

/"OO 

ll^ll 2 = (v\v) = / dufj*(w)fj(uj) = 1. 
Jo 

Using the continuous-mode bosonic commutation relations 

[&(<*>),&*(<*/)] =S(u-u/) = \d(w),$(u') 
(the other commutators being zero), it follows that 

[aM.otfo]] =l=\d[fj],cP[fi] 



(1) 
(2) 

(3) 
(4) 

(5) 
(6) 



Now let us assume that the pulses are initially prepared in an entangled quantum state of the type of a two-mode 
squeezed vacuum state 



|* in ) = exp{g*a[r ? ] d[fj\ - H.c.} |0) = S (?) |0>, 



(7) 



with q — \q\e tiPq being the squeezing parameter. Most studies of continuous-variable systems in quantum information 
processing have been based on such states [p|j7| Jl^ , ^5| | . After having passed the devices, the output state of the pulses 
reads as 



|*out> - cxp {q*a'[rj\ d'[fj] - H.c.} |0) = S'{q)\0}, 



where 



(ko rj* {uj)a! {uj) , 



d'[fj\ = / dujf]*{uj)d'(u), 
Jo 

and the continuous- mode operators d'(u)) and d'(us) are given by Jj^j ] 

2 
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(8) 

(9) 
(10) 



i=l 



(11) 

(12) 



[see 0, equation (A5)]. In equations ([TT|), (12) the continuous-mode operators &i(ui) = a(ui) and di(u>) = d(ui) belong 
to the fields entering the first input ports of the two four-port devices, whereas 0,2(0)) and d2{uj) belong to the fields 
entering the second input ports, and the bosonic operators (ji(u)) and hi(u) describe excitations of the two devices. 
The quantities Tu(lu) = T(uj), Tn(ui) = T(w) and T 2 i(w) = R(u>), T 2 i(w) = R(u>) are respectively the transmission 
and reflection coefficients of the four-port devices with respect to the incoming fields at the first input ports. In 
what follows we assume that the second input ports of the devices are unused, that is, the fields there are in the 
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vacuum state, the corresponding variables together with the device variables being referred to as the environment £. 
Moreover, we assume that the devices are not excited. 

It is convenient to represent each of the operators a'[rj] and d'[fj] in equations (||), ( |l0| ) as a sum of two other 
independent bosonic operators, 

a'[ V ] = \\Tr,\\a[rf] + {l-\\Tr,f) 1 ^q ri , (13) 
£[fj] = \\ffj\\d[fj'] + (l-\\ffi\\ 2 y/ 2 p fj , (14) 

where rj'(u>) = T(u))r)(uj)/\\Tr]\\ and fj'(uo) — T(u>)fj(u))/\\Tfj\\. Note that the operators q v and fifj belong to the 
environment £ , 

2 

(1 - H^ll 2 ) 172 % = a 2 [R V } + J29i[FiiV}, (15) 

i=l 
2 

(1 - llfr/ll 2 ) 1 / 2 ^ = d 2 [m] ■ (16) 

i=l 

Given the output quantum state |\& ut) of the system, the (symmetric) characteristic function 

*out(a, P) = (expfaatfo'] + /M^'] - H.c.) ) (17) 
can be calculated. In [B| it is shown that 

$out (a, /?) = exp { - i (|a| 2 [l + (cosh2|g| - 1) ||Try|| 2 ] 
+ |/f [l + (cosh2|g| -l)||f^|| 2 ]) 

-\ (a(3e- ilp « sinh2| g | + a*t3*e ilp " sinh2|g|) ||Trj||||f^||}. (18) 

Note that <E> out (a,/3) is a function only of the moduli |T(w)| and |T(w)| of the transmission coefficients. For T(u>) = 
T(lo) — 1, the characteristic function of the incoming fields is recognized. 

III. PULSE CORRELATIONS 

A measure of the entanglement of two subsystems (A) and (B) of a composed system (AB) that is prepared in some 



mixed state g is the quantum relative entropy E 28 1 , the quantum analog of the classical Kullback-Leibler entropy, 



E(g) = min Tr [g (In g—lxia)], (19) 

aeS 

where S is the set of all separable quantum states the composed system can be prepared in. For pure states, g = 
|5 , )(^E'|, the entanglement measure (|19| ) reduces to the von Neumann entropy of one subsystem 

E(g) = S A = Tr^HgWlngW) = S B , = Tt^g (20) 



[X r (^) (Tr^ s '), trace with respect to the subsystem A (B)]. In the case where |^) is the two-mode squeezed vacuum 
state (Q) it follows that 

E{g) = 5 t h(n sq ) = (n sq + 1) In (n sq + 1) - n sq lnn sq , (21) 



with n sq = sinh 2 |g| being the mean number of photons in each mode. Hence, the entanglement is given by the von 
Neumann entropy of a thermal state of mean photon number n = n sq . Note that for large mean photon numbers the 
entanglement increases linearly with the squeezing parameter, E w lnn sq w A\q\. 

Unfortunately, there has been no explicit expression for calculating the entanglement measure (|l9| ) of mixed states, 
which are typically observed in noisy systems. Extensive numerical procedures would be indispensable in general, 
which dramatically grow up with increasing dimension of the Hilbert space p9|| . Therefore some other correlation 
measures have been introduced. Although they are not purely quantum correlation measures [as is the entanglement 
measure ([l9h], they may be very helpful to gain insight into the problem of degradation of quantum correlations. 
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A. Fidelity 



Fidelity can be regarded as being a measure of how close to each other are two (system) states in the correspond- 
ing Hilbert space. Fidelities have been used, e.g., to describe decoherence effects in the transmission of quantum 
information through noisy channels 1 30 3l[] and to characterize the quality of quantum teleportation S| . 

When {3i n = l^in) (^inl and g out = I'J'outX^outl are respectively the input and the output density operators of the 
overall system (composed of the two pulses and the environment), then the density operators of the incoming and 
outgoing fields are respectively 



and 



fout ~ ±L fout ) 

(Tr( £ \ trace with respect to the environment). Following c, the fidelity 



(22) 



(23) 



F e = Tr(gt d) ^ t d) ) (24) 
can be defined. In O it is shown that applying the quantum state transformation outlined in Section O leads to 



F F = 1 



1 - {r}\Tr}){fj\ffi) n. 



' S( | 



(25) 



Note that F e in equation ( |25| ) refers to the transmitted light. Replacing T(u>) and T{lo) with R(uj) and R(u>) 
respectively, the fidelity with respect to the reflected fields is obtained. From equation ( J25| ) it is seen that the 
fidelity sensitively depends on the spectral overlaps of the transmitted and the incoming pulses, and these overlaps 
are substantially determined by the dependence on frequency of the transmission coefficients T(uj) and T(u>) (cf. 
Fig. ^). It is worth noting that the fidelity decreases rapidly with increasing initial mean photon number, that is, 
with increasing initial squeezing and thus increasing initial entanglement. 



q\ = 0.1 



q\ = 0.5 



«| = 2 

-1 -0.5 0.5 1 

(ri\Trj)(ij\fij) 

FIG. 1. The fidelity F e , equation is shown as a function of (real) (r]\Tri)(fj\frj) for various values of the squeezing 

parameter \q\. 

To illustrate the effect, F e is shown in Fig. |l| as a function of the product of the overlaps (r)\Tr))(f}\Tf}) for different 
values of the squeezing parameter \q\. In the figure (r]\Tr])(fj\Tfj) is assumed to be real. If the phases of T(w) and 
T(lu) did not depend on uj, then real MTi^)(fi\Tf]) would correspond to the maximally attainable fidelity (which could 
be also called entanglement fidelity ||32|| ). 

In quantum teleportation a strongly entangled two-mode squeezed vacuum is desired. Since in this case the photon 
number must be large, the state should tend to a macroscopic (at least mesoscopic) state, and hence its nonclassical 
features can become extremely unstable. Clearly, in the case of discrete-variable systems the entanglement must not 
necessarily increase with the mean photon number. For example, in the case of a two-mode state of the type 
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I*in) = 



(|00} + A|nn)) = 



1 + A 



(a[rj\ d[fj\) 



10) 



the entanglement is 



E 



Sth(|A| 2 ) 
1 + |A| 2 ; 



(26) 



(27) 



whereas the mean photon number in one mode is n — n|A| 2 /(l + |A| 2 ). The entanglement attains its maximal value 
E = In 2 at |A| = 1, which corresponds to a Bell-type state. Nevertheless, when n becomes large the entanglement of 
the transmitted field decreases exponentially with n J29|. A similar behaviour is observed for the fidelity. Using the 
results in Ref. [129], the fidelity (with respect to the transmitted fields) is obtained to be 



F„ 



1 



(1 + |A| 2 ) 2 



f + \X\ 2 (r,\Tr 1 r(fj\ffjr 



+ |A| 2 (l-||rr / || 2 )"( 1 -M 2 )' 



It is seen that with increasing value of n the fidelity rapidly decreases to the minimal value, i.e., F e 



(28) 

0.5 for |A| = 1. 



B. Entropic correlation measures 

Correlation measures can be defined employing the von Neumann entropy. A very general correlation measure is 
the index of correlation |33| 

Ic = S a + S d - S ad , (29) 

where 

S ad = -Tr(<tf ln^ u f) (30) 
is the entropy of the two-pulse system, and 

S ro = -Tr(^ t ) lnp(™ ) ) (31) 
(to = a, d) are the entropies of the single-pulse systems, 

_ rryCd) ~(a,d) z(d) _ mfa) ~(a,d) 

fout ^ 11 ^out > ^out — 11 tfout ■ l JZ 7 

Note that I c is bounded from below by \S a — Sd\ and from above by S a + Sd- Further, I c is an upper bound of the 
entanglement, E < I c , because g^ t §5 g^ t is an element of the set of separable states S of the two-pulse system. 
Correlation measures that can be used to formulate criteria of nonclassical correlation are pO| , pT| 

le ^ = £*m — S a d (33) 

(in = a,d). Since the entropy of a classical system must not be less than the entropy of one of its subsystems, positive 
values of li m ^ indicate nonclassical correlation. Thus, positive values of Ie and/or 1^ may be regarded as indicating 
entanglement. 

As can be seen from equation (Il8t), the characteristic function of is of Gaussian type. The same is true for the 



characteristic function of pj°2 and g^ t , which follows from equation (|18|) for = and a = respectively. Hence, 
the entropies S m , equation (pl|), and [after diagonalizing the quadratic form in the exponent in equation (fl8|)1 the 
entropy S a d, equation (|30|), can be obtained analytically in the form of the entropy of thermal states: 

S a = S th {n a ), n a = ||T?7|| 2 n sq , (34) 
Sd = S t h(rid), n d = ||T77|| 2 n sq , (35) 

Sad = S t h(n a d), n a d = (1 - ||rr?|| 2 ||T^|| 2 )n sq . (36) 

Equations ( |34] ) - ( ^6| ) again reveal the typical dependence on the spectral overlaps of the transmitted and the incoming 
pulses, the overlaps being substantially determined by the frequency response of the devices. Note that only the 
absolute values of T(ui) and T(uj) appear. Obviously, the entropies of the reflected light are obtained by replacing 
T(u>) and T(w) with R (uo) and R (cu) respectively. 
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FIG. 2. The entropy .Sad, equation (pOj), and the cor- 
relation indices I c , equation p9|), and /i and 7e , equa- 
tion (fj3|), are shown as functions of ||T?7|| 2 = (Trj\Tri) for 
\q\ = 2 [f (u) = 1]. 
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FIG. 3. The normalized correlation indices I e (solid 
line) and 7~i d ' (dashed line), equation (^), are shown as 
functions of ||Tt?|| 2 = {Tr\\Trj) for various values of the 
squeezing parameter \q\ [T(u>) = 1]. The values \q\ = 0.1, 
\q\ = 2, and |g| = 3, respectively, correspond to the mean 
photon numbers n sq = 0.01, n sq = 13, and n sq = 100. 



Examples of the entropy S a d and the correlation indices I c and Ie as functions of HT^H 2 are shown in Fig. || for 
||f?7|| 2 = 1 (that is, only one channel is noisy). It is seen that for not too small values of ||T?7|| 2 both li a ^ and 1^ 
are positive and thus indicate entanglement. With decreasing value of ||Tr/|| they decrease in a similar way as I c . 
Whereas Ie referring to the noise channel becomes negative and thus attains classically allowed values, 1^ referring 
to the unperturbed channel remains always positive but becomes small. The results are in agreement with the Peres- 
Horodecki separability criterion for bipartite Gaussian quantum states Jig] , It tells us that in the low-temperature 
limit the two-pulse system under consideration remains inseparable for all values of ||Tr7|| 2 and ||T7)|| 2 (cf. p4||). In 
order to get insight into the influence of the initial mean photon number n sq — sinh 2 |g| on the degradation of the 

quantum correlation, we have plotted in Fig. || the dependence on ||Try|| 2 of the correlation indices ii^ and ii^ for 
different values of \q\. It is clearly seen that the larger h sq becomes, the faster /| a ' and 1^ decrease. 



C. Example: absorbing dielectric plate 

The values of (t]\Tt]), (fj\Tf)), and 1 1 ZZ""-^ 1 1 2 , ||T7)|| 2 depend on the chosen pulse forms and on the dependence on 
frequency of the transmission coefficients of the devices T(u>) and T(u). When the bandwidths of the pulses are 
sufficiently small, so that the variation of T(u>) and T(uj) within the pulse bandwidths can be disregarded, then T{uj) 
and T(lu) can be taken at the mid-frequencies ujf and Qf of the pulse, i.e., (j]\TijTj) ~ Tij{u>f) and [fl\Tijfj] « Tij(ujf). 
In this case, the results become independent of the pulse shape and solely reflect the effect of the devices. 

Let us again consider the case where one pulse passes through free space and assume that the other pulse passes 
through a dielectric plate whose complex (Lorentz-type) permittivity is given by 



1 — (uj/ujo) 2 — 2ijLo/ujQ 

Explicit expressions for the transmission coefficient T(ui) and the reflection coefficient R(oj) of such a device are given 
in Ref. p5fl . In Figs. [|(a) and (b) they are shown as functions of to for e s = 1.5, j/loq = 0.01, and plate thickness 
2c/loq. Typically, the transmission amplitude shows a dip near the medium resonance ljq, whereas the reflection is 
peaked there. In Figs. ||(c) and (d), the fidelities F e of the transmitted and reflected fields are shown as functions of 
u)f. Finally, Figs. |](e) - (h) show the dependence on w/ of the corresponding correlation indices Ie . In particular, 

it is seen that for transmission F e and Ie are strongly reduced near the medium resonance, ujf ~ uoq 1 whereas for 
reflection they are enhanced in that region. 
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FIG. 4. The fidelity, equation (B5h 



2 3 
frequency 

and the correlation indices 1^ and 7^ > equation (|3^) are shown as functions of the 
mid- frequency u>f of narrow-bandwidth pulses that are transmitted (left figures) and reflected (right figures) at a dielectric plate 
for various values of the squeezing parameter \q\ [T(lj) = 1]. For comparison, the dependence on frequency of the amplitude and 
the phase of the transmission coefficient T(io) and the reflection coefficient R(ui) of the dielectric plate [complex permittivity 
( p7[ ) with e s = 1.5, j/uJo = 0.01; thickness 2c/luq\ are shown. 



IV. CONCLUSIONS 



We have studied the problem ol quantum correlations in a system of two entangled optical pulses of arbitrary 
shape, propagating through dispersing and obsorbing four-port devices, which may serve as models of beam splitters, 
fibres, interferometers etc. Given the complex refractive-index profiles of the devices, which can be determined 
experimentally, the relevant parameters for the pulse propagation and the associated quantum-state transformation 
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can be calculated without further assumptions. The devices can be viewed as realizations of noisy channels for 
quantum information transmission. 

In the calculations we have assumed that the pulses are initially prepared in a two-mode squeezed vacuum, which is 
typically considered in quantum teleportation of continuous- variable systems. We have calculated various correlation 
measures and studied their dependence on the pulse shape and mid- frequency, the frequency response of the devices, 
and the initial mean photon number. In particular, the results suggest that nonclassical correlations rapidly decrease 
with increasing mean photon number, which may drastically limit the effectively realizable non-classical correlation 
in continuous-variable systems. As a realization of a four-port device, we have considered a dielectric plate of a 
Lorentz-type complex permittivity in more detail. 

Finally, it should be remembered that the fidelity and the correlation indices are not strict entanglement measures. 
Nevertheless, they may be helpful to find characteristic dependences and to estimate limits of entanglement transmis- 
sion. Whereas for low-dimensional discrete- variable systems such as qubits numerical methods can be used to calculate 
the entanglement degradation exactly, the exploding effort prevents one from applying them to continuous-variable 
systems such as strongly squeezed two-mode vacuum states. It has therefore been a great challenge to find explicit 
entanglement measures, at least for some classes of states such as Gaussian states. 
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APPENDIX A: INPUT-OUTPUT RELATIONS 



The operator input-output relations at a dispersive and absorbing four-port read p5| , |35| ] 

2 2 

t>j M = T i i ^ &i ^ + A ^ H 9i M . ( A1 ) 

i=l i=l 

where a,i(u>) and bj(u>) are respectively the input- and output (destruction) operators of the radiation, and <7,(w) are 
the bosonic (destruction) operators of the device excitations. The 2x2 matrices T(ui) (transformation matrix) and 
A(w) (absorption matrix), which are determined by the complex refractive-index profile of the device pEfl , satisfy the 
relation 

T(w)T+ (tu) + A(uj) A+ (w) = I, (A2) 

provided that the device is embedded in vacuum. Whereas the matrix T(w) describes the effects of transmission and 
reflection, the matrix A(u) results from the material absorption. 
Let 

Qin = Qm[a(u>),a\oj)] (A3) 

be the overall input density operator, which is an operator functional of dc(u>) and a (u>), with a(us) being a "four- 
vector" according to 

a(uj) = [ai(u), a 2 (uj), gi(u>), g 2 {u)] T . (A4) 
The operator input-output relation flAip is then equivalent to the quantum-state transformation 19 

g out = Q ia [A + (u)6c(u),A T (w)at i (u)]. (A5) 
The unitary 4x4 matrix A(u) can be expressed in terms of the 2x2 matrices T(w) and A(ui) as 

/ T( W ) A(v) \ 

V F (^) G M / 

where F{uj) = -S(w)C- 1 (w)T(w), G(lu) = C(w)S- 1 (w)A(w), and C(cu) = v /T(w)T+(w), S(w) = ^A(w)A+(w) are 
commuting positive Hermitian matrices with 

C 2 (w) + S 2 (w) = 1. (A7) 
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APPENDIX B: DERIVATION OF EQUATION (18) 



Combining equations (ph and ( |17[ ) yields 

*out(a, /3) = (exp (afiV] + - H - 



S 



n (q) exp (aa f [77'] + 0[rj'] - H.c.) S" (g) 





exp^aS't (g) tf[rf] S" (g) + /3S* (g) tftft'] S" (g) - H.c 
It is not difficult to prove, on applying the operator expansion theorem, that 

S* (q)a[r,'] S' (g) = 
= o[i/] - sinh|g|||rr ? || (||Tr)|| <^'] + (1 - ||f ^|| 2 ) 1/2 pt 

+ (cosh |g|-l) ||T77|| (UTtjII &[t/] + (1 - \\T V \\ 2 ) 1/2 %) 

and 

5't (q)d[fj'] S' (g) = 

= d[ifl - e*« si n h|g|||ffy|] (||Try|| a^'] + (1 _ \\Tr,\\*)V* q% 
+(cosh|g|-l) \\ffj\\ (\\ffj\\ d[fj'} + (1 - WffiWY^Pv 



(Bl) 



(B2) 



(B3) 



Substituting these expressions into equation (Bl) and performing a normal ordering procedure, after some straightfor- 
ward calculations we arrive at <& ou t(a,/3) as given in equation (|l8|). Note that for T(ui) — T(uj) = 1 the characteristic 
function $i n (a, 0) of the quantum state of the incoming fields is obtained. 



APPENDIX C: DERIVATION OF EQUATION (|25|) 

According to equations (fit), (||), and (p4|), we may write 

p _ Tr (M) / s(a,<*)i(a,d)\ 
r e — -Ll I J)j n ^ out I 

= Tr( Q ' d ' {£(g)|0 o , d )(0 O) d |5t(g)Tr^ [s"(g)|0)(G|,S't(g)] } 
= Tr{s' Hq)S(q)\O a , d )(0 a , O d |£t(g)S'(g)|0)(0|} 
= (0|5'' t (g)^(g)|0 a , d ){0 a , d |5'+fQ)^'(9)|0> . 
We introduce the series expansions 



%)|0a,0 d ) = ^^^ tanhn M e " 



■impq 



K,o d ), 



(Cl) 



(C2) 



n=0 



and find 



5'(g)| ) = V tanh™|g|e— ^ |0) 

cosh g ^ 11 to! 1 ' 

m m = 



00 . 

(0 Q , d |5t (g) 5' (g) |0) = —j- £ tanh 2 "|g| 
cosh |g| ^ I 

x (U a , 0d| : : |U) 



(C3) 



(C4) 
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Combining equations (CI) and (C4) and using the the commutation relations 

[afoUV]] = (v\Tv)/\\Tvl 

'dmMv']] =(fj\ffj)/\\ffj\\, 

we eventually arrive at 



(C5) 
(C6) 



F K = 



1 



cosh \q\ 
1 

cosh 4 |g| 



J2 [(v\T V )(v\fv)tanh 2 \q\ 



rn—0 



1 -WTrjmf f,)tanh \q\ 



= 1+ 1 - (r)\Tr))(fi\ffj) smh 2 \q 



(C7) 
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